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I. INTRODUCTORY MATERIAL 



A. Groups and Algebras 

1. Definitions 



The mathematical tool one uses to study symmetries of physical systems is the theory of groups and algebras. A 
group is a set G = {a, b, c, .... } on which a multiplication operation is defined with the properties: 

• If a & b are in G, a b is also in G. 

• There is an identity element e: e Q a — a e = a for any a in G. 

• For every a in G, there is an inverse element in G, called a -1 such that a a -1 = a -1 a = e. 

• For every a, b, and c in G we have (a b) c = a (b c). 

For an Abelian group this operation is commutative: a b — b a. A group is continuous if its elements are 
functions of one or more continuous variables. A group is called continuously connected if a continuos variation of its 
variables leads from one arbitrary element of the group to another. Such groups are Lie groups. Lie groups whose 
parameters range over closed intervals are called compact Lie groups. Two groups a, b, c, ... and a', b', c', ... arc 
called isomorphic if a bijective transformation between elements of both groups exits (a ■<-> a, , b O b' , ...) such that 
a0if> a' 6', etc. Finally if a group G\ is isomorphic to another group G2, whose elements are matrices, G2 is 
called to be a matrix representation of G\. 

2. O(N) and SO(N) 
Consider a column vector in the N-dimensional real space 

( ^ \ 

x = X2 with the norm x T x — x\ + x\ + ... + x 2 N . (1) 

V XN J 

O(N) is the group of transformations, x — > x' — Ux, which leave this norm invariant: 

x lT x' = x T U T Ux = x T x U T U = 1 (detW = ±1). (2) 

Hence O(N) is the group isomorphic to the group of N x N, real, orthogonal matrices. If one chooses matrices with 
detU — +1, then one gets the SO(N) subgroup. 
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3. Unitary groups 

U{N) is the group of transformations which leave the norm x*x — x\x\ + x\xi + ■•• + x* n xm in the N-dimensional 
complex space invariant, i.e. it is the group isomorphic to that of the N x N complex, unitary matrices : WlA = 1. 
Clearly one has detU = ±1. If we take only those matrices with detU = +1, then we get the SU(N) subgroup. 
SU(2), for example is the group composed of 2 x 2 complex matrices of the form 

u= {ii ~tl ) with i^i 2 +i^i 2 =i (3) 

This is the 2-dimensional (j = 1/2) representation familiar from quantum mechanical spin. There are also higher- 
dimensional representations. For example, the three -dimensional (j = 1) representation is 

p (i) = V2~rii>2 \^i\ 2 -m 2 -V2^r 2 (4) 

Note the one-to-one correspondence between these representations: for each parameter set ipi and tp2, there is one 
unique 2x2 and one unique 3x3 matrix. 



B. Lie Algebras and Lie Groups 

Consider a Lie group whose elements U(9i) are parameterized by variables 6i such that U(B\ = 0, Qi = 0, 6*3 = 0, ....) 
is the identity element (X) of the group and that 



(5) 



0=0 



(6) 



The generators of the Lie Algebra are defined to be 

Bt = l w 

u ° l all 0=0 

where the infinitesimal change from the identity in the i-direction is 

U{S9i)=l-iBid0i. (7) 
Hence the finite change from the identity in the i-direction is 

U{6i) = {1- iBi59i){l - iB l 56 l )...{l - ^B ^ <W ^ )• (8) 
Taking the limit as N —> 00 we get U(6i) — exp(—iBi9i). Hence quantities of the form 

U(8 1 ,0 2 ,...) = e~^ B ^ (9) 
form a group if the following condition is satisfied: 

e-'E^-^^E,^ = e -iJ2i B i e i( e ^\ (10) 



The question is what restrictions need to be imposed on the quantities Bi to satisfy the condition in Eq. (|10p. The 
answer is given by the Baker-Campbell-Hausdorf Lemma, which states 

e A e B = exp(^A + B + i[A,B] + l[[A,B],A]] + l[[A,B],B]] (11) 

— —[[[A, B], B], A] + more nested commutators^ (12) 

Hence if a set of operators close under commutation relations {this is the definition of a Lie algebra), then they 
generate a Lie group: i.e. if [Bj, Bj] ~ B^. then 

Elements of \ _ / / continuous \ / Elements of 
Lie group J 6 I I parameters J I Lie algebra 



In many problems of physics one needs the evolution operator, U: 

ifi^- = HU with lim U(t) = 1. (13) 

at t->-oo 

It follows form the discussion above that if the Hamiltonian H is a sum of the elements of a Lie algebra, then the 
evolution operator U is an element of the corresponding Lie group. 

II. REALIZATIONS OF THE LIE ALGEBRAS 
A. Matrix Realizations 

For the SU(2) algebra, [Jj, Jfc] = icijk^k, the lowest dimensional representation (also called fundamental represen- 
tation) is provided by the Pauli matrices: 

For SU(N), N > 3, the N x N matrices, realizing the lowest (N-dimensional) representation are written as 

where Aj are usually referred to as Gell-Mann matrices. 

1. Fock Space Realization - Bosons 
For SU(N), introduce N boson creation and annihilation operators: 

[bi,bt]=5 tf) [b j ,b J -]=0 = [bt,bt], i,j = l,-..,N. (15) 
It is straightforward to show that the operators 

B » b J -,a = l,...,JV 2 or(7V-l) 2 (16) 

i,j=l ^ ' ij 

satisfy the same commutation relations as (A /2). This is called a change of basis of the algebra: 

B a Tij = bjbj with [T^ , T kn ] = S jk T in - 5 in T kj . (17) 

2. Fock Space Realization - Fermions 
For SU(N), one can also use N fermion creation and annihilation operators: 

{a Q , a^} = 5 a p, {a Q , a^} = = {a^,a^}, a, = 1, • • • , N. (18) 

The operators 

Q af3 = aj, a fj (19) 
also satisfy the same commutation relations as T,j . We then have two different representations of the SU (N) algebra: 

Tij : completely symmetric representation, 

Q Q/ 3 : completely antisymmetric representation. 
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B. Invariants and Labeling 

1. Casimir Operators 



A Casimir operator, C, is an operator which commutes with all the elements of the algebra: [C,B a ] = 0. Schur's 
lemma states that C oc I . For U(N), there are N independent Casimir operators: 

Ci = Tn, 

i 

C 2 = y^TjjTjj, 

i,i 

C3 = TjjTjfeTfcj, 

i,3,k 



2. Labeling of States 

To write down the states associated with a given algebra we first need to find all the subalgebras included in this 
algebra: 



Ad AiD A 2 d 



Then the state can be written as 



(state) = |ai,a 2 ,..;/3i,^ 2 ,-;7i)72,- : •••) 

a%, a-2 : eigenvalues of Casimirs of A 
/?i , /?2 : eigenvalues of Casimirs of Ai 
71 , 72 : eigenvalues of Casimirs of A2 ■ 

For example, for the SU(2) algebra the familiar projection on the third component of the angular momentum yields: 

SU(2)DU(l)^\j,m). 

C. Example: Boson Fock-Space Realization for SU(2) 

Consider the boson realization of the SU(2) algebra given in Eq. (fT6|) : 

J+ = Ji+iJ 2 =blb 2 , J- = (J+) 1 =b 2 b 1; (20) 



J 3 = -(b t 1 b 1 -btb 2 ). (21) 

To find out what values of the representation labels are permitted one should write down the Casimir operator in 
terms of particle creation and annihilation operators: 

1 N /N \ 

c 2 = J§ + ^(J+J- + J-J+) = y [j + ij = i(i + i), (22) 

where 

N = b|b!+b 2 b2. (23) 
Hence fixing the number of particles fixes the label j. 
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III. QUASIPARTICLES AND BOGOLIUBOV TRANSFORMATIONS 

Consider one fermion states 

\a) =a+|0), a Q |0) =0 

and the associated time-reversed states: 

|-a)=aL a |0), a_ a |0> =0 
Note that and al Q are distinct and they anticommute. Quasi-particle operators are defined as 

A Q = u a a a - v a a}_ a , (24) 

A- a = U a a- a — VaB^. (25) 

Ajj and AL Q are obtained by Hermitian conjugation. Imposing the condition that the quasi-particles satisfy the 
canonical anticommutation relations yields |u a | 2 + |t» Q | 2 = 1. We can rewrite the connection between quasi-particle 
and particle operators as 



Note that this is an SU(2) group transformation. Corresponding SU(2) algebra is called quasi-spin algebra [T|: 

S«=ala t _ Q , S«=(S«)t (27) 

S^^a^ + al^-l). (28) 
There are as many commuting SU(2) algebras as the possible values of a: 

[S£,S*] = 2S£<5 Q/3 , (29) 

[S£,s£] = ±S£<T' 3 . (30) 

Clearly the realization in Eqs. (f27| and (f28| is different than the realization given in Eq. (fl9l) . Indeed with N 
different fermion creation-annihilation operators 

a Q ,a^, a =!,■■■ N 
one can form the generators of the SO(2N) algebra: 

a a a a'' a aa Q ', a^a Q ' 

SU(N)subalgebra 

For N = 2, one gets 5*0(4) ~ SU(2) x SU(2). These two SU(2) algebras commute with one other; one of them is the 
algebra given in Eq. (|T5)) and the other one is the quasi-spin algebra. 

The transformation in Eq. (|26[) . called Bogoliubov transformation, can be written as an operator transformation 
under the quasi-spin SU(2) group: 

A a =TZa a TZ\ (31) 

where the group rotation is 

n = e - l Q s o ° e ,»s^ e io g (i+| 2Q | 2 )s o ° e ~,:s° (32) 

with 

z a = , e ^ = : . (33) 

U a \u\ a 
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Note that ground states (vacua) for particles and quasi-particles are different: 

|0) : Particle Vacuum 

\z) : Quasi — particle Vacuum 

with 

l*>=ft|0>. (34) 

and 

A a \z) = na^n^) = o. (35) 

IV. DYNAMICAL SYMMETRIES 

Consider a chain of algebras (or associated groups): 

G\ D G*2 D ■ ■ ■ D G n 

If a given Hamiltonian can be written in terms of only the Casimir operators of the algebras in this chain, then such 
a Hamiltonian is said to possess a dynamical symmetry: 



H = Y J [aiCi{G i ) + p i C 2 {G i )] 



i=i 



Obviously all these Casimir operators commute with each other, making the task of calculating energy eigenvalues 
straightforward. 

One of the more commonly used collective models in nuclear physics is the Interacting Boson Model 2]. In this 
model low-lying states of medium-heavy nuclei are obtained as states generated by six interacting bosons, one with 
angular momentum L = and five with angular momentum L = 2. It then follows from Eq. (|16p that bilinear 
products of the associated creation and annihilation operators form an SU(6) algebra and, if one limits the terms in 
the Hamiltonian to include at most two-body interactions, for certain values of the interaction strength one gets three 
dynamical symmetry chains: 

• Vibrational Nuclei: SU(6) D SU(5) D 5*0(5) D 50(3) 0, 

• Rotational Nuclei: SU(6) D SU(3) D 50(3) 0, 

• 7-Unstable Nuclei: SU(6) D 50(6) D 50(5) D 50(3) @. 

This model is covered in depth by other lecturers at this summer school @. 



A. Supersymmetry and Superalgebras 

1. Contrasting Symmetry and Supersymmetry 

As we have seen in the previous sections ordinary symmetries either transform bosons into bosons or fermions 
into fermions. Natural mathematical tools to explore them are the Lie groups and associated Lie algebras. We will 
designate a generic element of the Lie algebra as bosonic, Gb, symbolically: 

[Gb, Gb] = Gb 

Supersymmetries, on the other hand, transform bosons into bosons, fermions into fermions, AND bosons into 
fermions and vice versa. Tools to explore them are superalgebras and supergroups. A superalgebra is a set with two 
kinds of elements, Gb and Gf ■ It closes under commutation and anticommutation relations in the following manner: 

[Gb,Gb] = Gb, 
[Gb,Gf] = Gf, 
{Gf,Gf} — Gb- 

A simple example of a superalgebra is given in the next section. 
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2. Case Study: Simplest Superalgebra 
Consider three dimensional harmonic oscillator creation and annihilation operators ([6i,&J] = Sij) and define 



\i=l / i=l 



which lead to the commutation relations 

[K , K±] = ±A±, [K+, A_] = -2K . (37) 

This is the SU(1,1) algebra. It is non-compact. Recall that Casimir operators obtained by multiplying one, two, 
three elements of the algebra are called linear, quadratic, cubic Casimir operators. For SU(1,1) the quadratic Casimir 
operator is 

C 2 = Kl - X - (A+A_ + K_K + ) . (38) 

We next introduce spin (fermionic) degrees of freedom in addition to the bosonic (harmonic oscillator) ones and 
define 

i i 

One can show that the following commutation and anticommutation relations hold: 

[K , F±] = ± l -F ± , 

[K+, F+] = 0=[K_, F_], 
[K ± , F T ] = T F ± , 
{F±,F ± }=K ± , {F+,F_}=K . 

Along with the commutation relations given in Eq. p7[) . these are the commutation relations of an orthosymplectic 
superalgebra, Osp(l/2), which is also non-compact. Hence the operators K +1 K-, Kq, F + , and A_ are the generators 
the Osp(l/2) superalgebra. We have Osp(l/2) D SU(l,l). 
The Casimir operators of Osp(l/2) are given by 

C 3 (Osp(l/2)) = i(L+^) 2 = ij 2 , (40) 



4 V 2 

C 2 (SU(1,1)) = ^L 2 - A. (41) 



H=l( p 2 + r 2 )+\(a-L + ?) (42) 



Hence a Hamiltonian of the form 

11 = 

2 vr ' \ 2 

can be rewritten in terms of the Casimir operators of the group chain Osp(l/2) D SU(1, 1) D SO (2) 0]: 

H = 4AC 2 (Osp(l/2)) - 4AC 2 (SU(1, 1)) + 2A (43) 

This is an example of a dynamical supersymmetry. In a microscopic interpretation the bosons of the Interacting Boson 
Model are taken to be correlated pairs of nucleons [8j . For odd-even nuclei, the algebraic structure of the Interacting 
Boson Model can be extended to include the unpaired fermions. In such extensions dynamical supersymmetries 
naturally emerge Q. In this case unpaired fermions in Ji,j2,j3,''' orbitals can be placed in a fermionic algebra 
of SUF(J2i(^ji + !))• The resulting SU(6)b X SUF(J2 i (2j i + 1)) algebra is then embedded in the superalgebra 
SV (6 / SU F(%2i(2ji + !))• There are several experimental examples of such dynamical supersymmetries |10l - tl2j |. 
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B. Supersymmetric Quantum Mechanics and Its Applications in Nuclear Physics 



Consider two Hamiltonians 

Hi =G t G, H 2 = GG\ (44) 
where G is an arbitrary operator. The eigenvalues of these two Hamiltonians 

G^G\l,n) = EM\l,n) 
GGt|2,n) = E^\2,n) 

are the same: 

= = E n (45) 

and the eigenvectors are related: |2, n) = G [G^G] |1, n). (This works for all cases except when G|l, n) = 0, which 
should be the ground state energy of the positive-definite Hamiltonian Hi). 

The pair of the Hamiltonians in Eq. f|44[) define the supersymmetric quantum mechanics [131 ]. To see why this 
construction is called supersymmetry we define the operators 



Then it is easy to see that the "Hamiltonian" 



0\ (0 G 

G* ' W loo 



B={Q,Q*} = ( H j I) (46) 

is an element of a simple superalgebra along with with the operators Q and 

[H,Q]=0 = [H,Ql]. 

Very few realistic Hamiltonians can be cast in the form given in Eq. (|46p . (A couple examples are given below 
below). However, sup ersymmetric quantum mechanics can be a starting point for a semiclassical expansion of most 
Hamiltonians flil Il5|. 

The nuclear shell model is a mean-field theory where the single particle levels can be taken as those of a three- 
dimensional harmonic oscillator (hence labeled with SU(3) quantum numbers) for the lowest (A < 20) levels. For 
nuclei with more than 20 protons or neutrons, different parity orbitals mix. The Nilsson Hamiltonian of the spherical 
shell model is 

H = ujbjbi - 2/sL.S - fc^L 2 , (47) 

where the second term mixes opposite parity orbitals and the last term mocks up the deeper potential felt by the 
nucleons as L increases. 

Fits to data suggest /x s» 0.5, which lead to degeneracies in the single particle spectra. In the 50-82 shell (whose 
SU(3) label or the principal harmonic oscillator quantum number is N = 4), the S1/2 and d-3/2 orbitals and further 
d 5 / 2 and g 7 / 2 orbitals are almost degenerate. It is possible to give aphenomenological account of this degeneracy 
by introducing a second SU(3) algebra called the pseudo-SU(3) [16|, [l7j- Assuming that those orbitals belong to the 
N = 3 (with I = 1,3) representation of the latter SU(3) algebra one designates the quantum numbers of the SO(3) 
algebra included in this new SU(3) to be pseudo-orbital-angular momentum (£ = 1,3 in this case) and introduces 
a pseudo-spin (s = |). One can easily show that j = 1/2 and 3/2 orbitals (and also j = 5/2 and 7/2 orbitals) 
are degenerate if pseudo-orbital angular momentum and pseudo-spin coupling vanishes. It was later discovered that 
pseudo-spin symmetry has a relativistic origin (l8| . 

It was shown that two Hamiltonians written in the SU(3) and the pseudo-SU(3) bases are supersymmetric partners 
of each other [l9( • The operator that transforms these two bases into one another is 

U = G [&G] ~ V2 = V2F_ (K + [F + ,F_}y 1/2 
= fai) (b\bi-aiL^~ 12 

yielding 

F pseudo-SU(3) = U i/ SU(3) C/t 
= b\bi - 2k (2/x — 1) L • S - kfiL 2 + [1 - 2k(ji - 1)] . 
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V. INFINITE ALGEBRAS 

Sometimes the algebras associated with the symmetries of the physical problems have an infinite number of elements. 
One example is an algebra originally introduced by Gaudin in his study of spin Hamiltonians [20] : 

[J+(A) ; J-(m)]=2 J0 ^- J0M , (48) 

A — /i 

[J°(\),J±M]=± J±{X) - J±M , (49) 
A — (j, 

[J°(X),J°(n)] = [J ± (X),J ± (fj,)}=0. (50) 

In the above equations A is an arbitrary complex parameter. To see the relevance of the Gaudin algebra to nuclear 
physics, let us rewrite the quasispin algebra of Eqs. (|2"T|) and (|2"5|) using fcrmion operators of the spherical shell model: 



§ t = E(- 1 ) ( ^ m)a ]™«]— ( 51 ) 

m>0 

SJ- = ]r(-l) (j - m) aj -™%™, (52) 



m>0 



m>0 

Note that these are mutually commuting SU(2) algebras: 



§j — — E ( a ] m a j m + a ] - m a j -m 1) ) ■ (53) 



A possible realization of the Gaudin algebra can be given in terms of the elements of the quasi-spin algebra (see e.g. 
Ref. (HH): 

N OQ N q± 



, £i - A f-( ei - A ' 

1=1 2—1 

where are arbitrary constants. The operator 

H(X) = J°(\)J°(\) + ij+(A)J"(A) + ±J-(\)J+(X) (55) 
is not the Casimir operator of the Gaudin algebra, but a conserved charge: 

[H(X),H(»)}=0, A^/x- (56) 

Lowest weight vector is chosen to satisfy 

J-(A)|0)=0, and J°(A)|0) = W(A)|0), (57) 

so that 

ff(A)|0) = [W(X) 2 -W'(X)] |0), (58) 

where prime denotes derivative with respect to A. It is also possible to write bosonic representations of Gaudin-likc 
algebras (22| . 

To find other eigenstates of the operator in Eq. ((55]) we consider the state |£) = J + (£)|0) for an arbitrary complex 
number £. One gets 

[H(X), J+(0] = ^ ( J + (A) J°(0 - J + (0 J°(A)) . (59) 
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Hence, if W(£) = 0, then J + (£)|0) is an eigenstate of H(X) with the eigenvalue 

£a(A) = [W{Xf - W{X)] - 2j&. (60) 

Gaudin showed that this procedure can be generalized. Indeed a state of the form 

.-.,U>= J + (^i)J + (^) ■ ■ • J+(6)|0 > (61) 
is an eigenvector of H(X) if the numbers ■ ■ ■ ,£n G C satisfy the so-called Bethe Ansatz equations: 

" 1 

^(U^r-T for a=l,2,...,n. (62) 

Corresponding eigenvalue is 

E n (X) = [W(\f - W'(X)] 2 ± W W ~Y&\ (63) 

To make the connection to the spin Hamiltonians Gaudin studied one considers the limit 

nm(A-e fe )ff(A)=7^ = -2V^^. (64) 

Since the conserved charges commute for different values of the parameter, Eq. (|64p implies that H(X) and 7£fc can 
be simultaneously diagonalized: 

[H{X),H(^)] =0 =► [H(\),K k ] = (65) 
[7^,ft fe ]=0 (66) 



One can also show that 



and 



$^fti=0, (67) 



^e i 7e i = -2^S i -S J -. (68) 

i i^j 

Eqs. fM|l and (f6"8")l are the spin Hamiltonians considered by Gaudin. 

Richardson considered solutions of the pairing Hamiltonian using a different technique (23| . Below we derive his 
results using Gaudin's method. Note that the Gaudin Algebra of Eqs. (|48[). (|49]) . and ([50]) can be satisfied not only 
by the operators J (A), but also by the operators J (A) + c for a constant c. In this case 

H(X) = J(A) • J(A) H(X) + 2c • J(A) + c 2 (69) 

which has the same eigenstates. To exploit this fact we introduce new operators that we name Richardson operators: 

lim (A - e fc ) (H(X) + 2c • S) = R k = -2c ■ S k - 2 V Sfc ' Sj " . (70) 
It is straightforward to show that 

[fr(A)+2c-S,iZ*] = 0, [Rj,R k }=0. (71) 

One can also prove the identities 

Ri = -2c • ]T S fc (72) 

i k 
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and 

J2 ^ = - 2 E e ' c ■ s< - 2 E s ' • s j- ( 73 ) 

i i i^j 

Using the above results one is then ready to write down the eigenvalues of the pairing Hamiltonian given by 

H = E e i°5 m*i m-\G\J2 S+S- . (74) 

jm jj' 

Indeed choosing the constant in Eq. (|70| as 

c= (0,0,-l/2|G|) 

one can show that the solvable Hamiltonian of Eq. (j6"9")l can be written as 

ifi = £ ^ + i g i 2 (E - i g i E ^ + • • • > ( 75 ) 

2 2 Z 

which is the pairing Hamiltonian of Eq. (ITU) up to a constant. 

VI. PHYSICS APPLICATIONS 

A. Pairing problem in Nuclear Physics 

Pairing plays a very important role in nuclear physics. In previous sections we discussed how the quadrupole 
collectivity of medium-heavy nuclei can be represented by nucleon pairs coupled to angular momenta zero and two. 
Over the years considerable attention was paid to exactly solvable pairing Hamiltonians with one- and two-body 
interactions. The pairing interaction was first presented by Racah in LS-coupling scheme [2~4j and was generalized 
to the j j-coupling scheme [25j | . Here we confine ourselves to the s-wave pairing case as represented by the quasispin 
algebra of Eq. (f5T|) . (f52"j) and (|55|) . Note that solvable models with both monopole and quadrupole pairing also exist 
j26j . Exactly solvable cases so far studied for the monopole pairing case include 

• The exact quasi-spin limit [l|: 

H = -\G\J2S+S-. (76) 

jj' 

• Richardson's solution, discussed above, for the case when the single particle energies are added to the Hamilto- 
nian in Eq.® [H 

" Evw"» f; E- s :/-v ( 7? ) 

jm jj' 

• Gaudin's model [2(|, which is closely related to the Richardson's solution. 

• The limit with separable pairing in which the energy levels are degenerate (the one-body term becomes a constant 
for a given number of pairs) [27H29| : 

H=-\G\J2^c r S+S7,. (78) 

jj' 

• Most general separable case with two orbitals (30| . 
Introducing the operators 

S+{x) = Y- — f4-£+ S-{x)=Y- — t*-i5-&". (79) 

V ' 1 - \Cj\ 2 X 3 W ^1- \Cj\ 2 X 3 
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one can show that the state [27|, |28| 

5 + (0)^(zW)...5 + (zWi)|0> (80) 
is an eigenstate of the Hamiltonian in Eq. (|78|) with energy 



j k=l z k 



if the following Bethe ansatz equations are satisfied: 

N-l 



E — , 9 ,.w = — + £ JiV) jiv) m=l,2,...iV-l. (82) 

Similarly 



(iV) (iv) 1 (W) _ (AT) 

3 Zm Zm k=l(k^m) Zm Z k 



S + (x[ N) )S + (xi N) )...S+( X ^)\0) (83) 
is an eigenstate with zero energy if the following Bethe ansatz equations are satisfied: 



-0-/2 N 1 

E 1/|e , 2 _^v) = E (AT) _ (AT) foreve ^ ^ = 1,2,..., iV- (84) 



The states in Eqs. (f80|) and (f83|) are eigenstates of the Hamiltonian in Eq. (|78[) if available single-particle levels are 
at most half full. One can show that, if the single-particle levels are more than half full, the state 

.S-(*W)|0> (85) 

is an eigenstate with the same energy as in Eq. (|81j) if the Bethe ansatz equations given in Eq. (|82|) are satisfied [28[ . 
In Eq. ([83)) |0) designates the state where all single-particle levels are completely filled. 

It turns out that one can find an exact solution for the case where there are only two single-particle levels j30j . i.e. 
consider the Hamiltonian 

j§ = E M° - E § ? + E ^ ( 86 ) 

3 33' 3 

where 6j and Cj's are dimensionless and the sums are performed over only two single-particle states. In the equation 
above we added a constant term for convenience where fij = j + -| is the maximum number of pairs that can occupy 
the level j. The eigenstates of the Hamiltonian in Eq. (f86|) can be written using the step operators: 

j 2£ 3 ~ \ C 3\ X 

as 

J + (x 1 )J+(x 2 )...J+(x N )\0). (88) 

Defining the auxiliary quantities 

P = 2 £jl — £j2 S = o £j2 ^ Cjl 1 2 " £jl ^ Cj2 ^ (89) 

\ c ji I 2 — I c i2 r l c ii I 2 — I c j2 r 



one obtains the energy eigenvalues as 

^ = -E^ (90) 



n=l 



In the above equations, the parameters Xk are to be found by solving the Bethe ansatz equations 

Slj\Cj\ 2 _ (3 | y, 2 ( 91 ) 

j 2e 3 - | Cj - \ 2 x k /3 - x k n J^ k) x n - x k 
A generalization of this approach to include three orbitals is still an open problem. 
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B. Subbarrier Fusion and Group Transformations 



In some applications one needs to use not only the algebra but also the entire group transformation. One such 
example is eikonal scattering from complex systems with dynamical symmetries [311 ] . Another example is provided 
by fusion reactions below the Coulomb barrier (32j . For fusion reactions near and below the Coulomb barrier the 
experimental observables are the cross section 



and the average angular momenta 



<w = 



The partial-wave cross sections in these equations are given by 



irft 2 
2^E 



(2l + l)T t (E), 



(92) 



(93) 



(94) 



where T[{E) is the quantum-mechanical transmission probability through the potential barrier and /j, is the reduced 
mass of the projectile and target system. The fusing system can be described by the Hamiltonian 



H = H k + V (r) + H (0 + H int {r,0 



(95) 



with the kinetic energy 



2[i 



(96) 



where r is the relative coordinate of the colliding nuclei and £ represents any internal degrees of freedom of the target 
or the projectile. In this equation Vo(r) is the bare potential and the term -ffo(£) represents the internal structure of 
the target or the projectile nucleus. 

The propagator to go from an initial state characterized by relative radial coordinate (the magnitude of r) ri and 
internal quantum numbers rii to a final state characterized by the radial position rf and the internal quantum numbers 
rif may be written as a path integral: 

K{r f ,n f , T-r i ,n h Q) = Jv [r{t)\ e^ r ' r V n/Tii (r(i),T), (97) 

where 5(r, T) is the action for the translational motion and Wn fni is the propagator for the internal system along a 
given path, [r(t)], of the translational motion: 



W njni (r,T) = {n f U in Mt),T) 



n, 



U- ln t satisfies the differential equation 



with the condition 



ih- 



dU ir 



Of 



[H + Hi n t] Uu 



(98) 



(99) 



U int (t = 0) = 1. 

In the limit when the initial and final states are far away from the barrier, the transition amplitude is given by the 
5-matrix element, which can be expressed in terms of the propagator as |32i | 



Si. 



lim 



P)P/ 



exp 



-(Pfrf-Piu) 



dTe +lET / h K{r f ,nf,T ]n , ni ,{)), 



(100) 
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where pi and p/ are the classical momenta associated with n and r/. In heavy ion fusion we are interested in the 
transition probability in which the internal system emerges in any final state. For the £th partial wave, this is 



which takes the form 



T e (E)=J2\Sn f , ni (E)\ 2 , 



(101) 



Tt(E) = lim 



PiPf 



cffexp 



V[r{t)] / X>[f(i)]exp 



—ET 
h 



Texp 



% -ET 

h 



~(S(r,T)-S(r,T)) 



Pm- 



(102) 



Here we have assumed that the energy dissipated to the internal system is small compared to the total energy and 
taken p/ outside the sum over final states. We identified the two-time influence functional as 



p M (f(t), T- r(t),T) = ]T W: ftni (f(t); T, 0)W n/ , nj (r(t);T, 0). 



(103) 



Using the completeness of final states, we can simplify this expression to write 

p M (r(t),T;r(t),T) = U uUr (t), f )U iat (r(t), T) 



(104) 



Eq. (|104[) shows the utility of the influence functional method when the internal system has symmetry properties. 
If the Hamiltonian in Eq. (|95l) has a dynamical or spectrum generating symmetry, i.e., if it can be written in terms 
of the Casimir operators and generators of a given Lie algebra, then the solution of Eq. (|M| is an element of the 
corresponding Lie group [33| . Consequently the two time influence functional of Eq. (I104[) is simply a diagonal group 
matrix clement for the lowest-weight state and it can be evaluated using standard group-theoretical methods. This 
is why the path integral method is very convenient when the internal structure is represented by an algebraic model 
such as the Interacting Boson Model. Using this approach it is possible to do systematic studies of subbarrier fusion 
cross sections [34[ as well as other observables [35j not only for nuclei that are described by the dynamical symmetry 
limits, but also for transitional nuclei. 



C. Neutrinos and their Symmetries 



The Standard Model does not contain neutrino masses. However, a neutrino mass term can be introduced as an 
effective interaction. Symmetries, in particular weak isospin invariance, define the Standard Model. In the neutrino 
sector this symmetry is SU(2)w x U(l). In the Standard Model, the left-handed and the right-handed components 
of the neutrino are treated differently: vl sits in an weak-isospin doublet (Iw =1/2) together with the left-handed 
component of the associated charged lepton, whereas vr is an weak-isospin singlet (lw=0)- A mass term connects left- 
and right-handed components. The usual Dirac mass term is L = = m(^> i^S r + ^r^l). But such a neutrino 

mass term breaks the weak-isospin symmetry, hence it is not permitted in the Standard Model. The right-handed 
component of the neutrino carries no weak isospin quantum numbers. As we will see below this feature permits 
Majorana neutrino mass in the Standard Model if one only uses right-handed neutrinos. 



1. Dimensional Counting 



Lagrangian, L, has dimensions of energy (or mass). In field theory one usually employs the Lagrangian density/!: 
L = J d?xC{x). Lagrangian density, £, has dimensions of energy / volume or M A . Usually one uses the terms 
Lagrangian and Lagrangian density interchangeably when the meaning is clear from the context, Defining a scaling 
dimension for x, [x] to be -1 then we see that the scaling dimension of momentum (or mass) should be [m] = +1 (recall 
that (p.x/K) is dimensionless and we take [7i]=0). Clearly one has [C] = 4. This should be true for any Lagrangian 
density of any theory. Considering the mass term for fermions, C m — mM/^ we conclude that \^>^] = 3 or [^>] = 3/2. 
In the Standard Model the Higgs field vacuum expectation value gives the particle mass: C = H^ffy, hence [H] = 1. 
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2. Effective Field Theories 

A Lagrangian describing a particular field should be consistent with the symmetries of this field, i.e. invariant 
under rotations, translations, Lorentz transformations, etc. (A combination of these symmetries describe the Poincare 
group). An example is provided by the Lagrangian of quantum electrodynamics. The two Lorentz invariants one can 
write down in terms of electric and magnetic fields are E 2 — B 2 and E • B. Under time reversal transformations the 
electric field does not change sign: 

E ->• E, 

but the magnetic field does: 

B B. 

Hence one of the Lorentz invariants (E • B) is not an invariant under time-reversal transformations. Requiring the 
time-reversal to be a good symmetry one writes the photon part of the Q.E.D. Lagrangian in terms of only the other 
invariant: 

£ 7 = i(E 2 -B 2 ). (105) 

The scaling dimension of the Lagrangian in Eq. (|105l) is, of course, four. In this problem there is a clear separation 
of energy scales. If the energy in the electromagnetic field is significantly below twice the mass of the lightest charged 
particle (electron), then there will be no energy loss to pair production. Effective field theories provide a framework 
to appraise the impact of the physics that takes place at higher energy scales on processes that occur at much lower 
energies. How do we take into account the physics at higher energy scales, e.g. the effect of the existence of charged 
particles on photons? (Essentially we are asking to integrate the charged particles out of the path integral for the full 
Q.E.D.). This effect can be represented by adding additional terms to introduce an effective Lagrangian: 

C -> ^effective = C + SC. (106) 

The additional Lagrangian should still be consistent with the symmetries of the system. Let us again use Q.E.D. 
as an example. Since we want the additional term to be Lorentz invariant, clearly it has to involve higher powers 
of Lorentz invariants. Since each Lorentz invariant for the electromagnetic field has a scaling dimension of four, the 
lowest dimensional (eight in this case) correction is 

a(E 2 -B 2 ) 2 + &(E-B) 2 , (107) 

where a and b are yet undetermined, dimensionless constants. Note that the square of E ■ B is time-reversal invariant 
even though E • B is not. However, the expression in Eq. (|107|) is still not a proper Lagrangian density since it does 
not have scaling dimension four. To make it four dimensional we need to divide it by some energy scale to the fourth 
power: 



a(E 2 -B 2 ) 2 + 6(E-B) 2 . (108) 



This is how far we can go with the effective field theory tools. However, an educated guess would suggest that the 
energy scale A should be proportional to the mass of the lightest charged particle in the leading order, A ~ m e . 
Directly integrating out the charged particle degree of freedom in the full Q.E.D. Lagrangian gives the numerical 
values of the dimensionless constants: 

C = i(E 2 - B 2 ) + -^-j [(E 2 - B 2 ) 2 + 7(E • B) 2 ] . (109) 
2 45m| 

This result is know as the Euler-Heisenberg Lagrangian in the literature [3(| . 

3. Neutrino Mass 



Even though the Standard Model does not include neutrino mass, it is possible to write effective Lagrangians for the 
neutrino mass in terms of Standard Model fields. Such a Lagrangian should preserve the SU(2)w x U(X) symmetry. 
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Recalling that 1^ — 1/2 for the vj, and —1/2 for iTsMj we can write a dimension-five operator describing neutrino 
mass using the Standard Model degrees of freedom: 

C = ^HsmHsm^vlp, (110) 

where z/£ Q is the charge-conjugate neutrino field and a and j3 are flavor labels. From the constants of Eq. (jllOl) one 
gets the usual neutrino mixing matrix 

^^ m (^~"W T (111) 

A y ' 

Clearly the neutrino mass term in Eq. (|110[) is not renormalizable. It is the only dimension-five operator one can 
write using the Standard Model degrees of freedom: In a sense the neutrino mass is the most accessible new physics 
beyond the Standard Model. 

A mass term of the type given in Eq. (|110j) is different than the usual charged-particle mass term in the Dirac 
equation and it is known as the Majorana mass term [37j . Such a mass term is permitted by the weak-isospin invariance 
of the Standard Model, but it violates lepton number conservation since it implies that neutrinos are their antiparticles. 
To gain a better insight into the nature of the Majorana mass term it is useful to consider transformations between 
particles and antiparticles. The particle- antiparticle symmetry, realized via the transformation 

*^a* + 6 75 * c , |a| 2 + |6| 2 = l (112) 

is usually referred to as Pauli-Gursey transformation (38[. It is easy to see that, under such a transformation, a Dirac 
mass term would transform into a mixture of Dirac and Majorana mass terms. The operators 

1 f J3, 



£>+ = -/ (PxVlVr, (113) 
L + = i / d 3 x , R + = \i d 3 x , (115) 



2 J 2_ 
their complex conjugates, and the operators 

Lo = \J d 3 x - , R = i J d 3 x - (116) 

form an SO(5) algebra [H]. 1 

The operators A +1 A_ and y4 = Rq — L form an SU(2) subalgebra that generates the Pauli-Giirsey transformation 
(SU(2)pq). The most general neutrino mass Hamiltonian sits in the SO(5)/SU(2)l x SU(2)r x U(1)l +r coset 
and can be diagonalized by a SU(2)pg rotation. This diagonalization is referred to as the see-saw mechanism in the 
literature (ilj . 

4- Neutrino Many-Body Theory 

Understanding neutrino propagation at the center of a core-collapse supernova requires a careful treatment of 
features like neutrino-neutrino scattering [42|, |43[ and antineutrino flavor transformations [3]. The neutrino self- 
interactions could especially impact the r- process nucleosynthesis taking place in core-collapse supernovae [45[ . There 
is an extensive literature on this subject, a good starting point is several recent surveys [46l. l47| . 

For simplicity, let us consider only two flavors of neutrinos: electron neutrino, v e , and another flavor, v x . Introducing 
the creation and annihilation operators for one neutrino with three momentum p, we can write down the generators 
of an SU(2) algebra 0: 

J+(P) = 4(P)«e(p), J-(P) = 4(P) a a(p), 

Mp) = \ (4(pK(p) - 4(pK(p)) ■ (117) 



1 Note the similarity between the Majorana mass term in Eq. I|110|l and the pairing interaction described by the quasispin algebra operators 
in Eqs. (15 1 I t and (1 52 I t ). Indeed the presence of an SO(5) algebraic structure is a general feature of particular pairing interactions |4(J . 
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Note that the integrals of these operators over all possible values of momenta also generate a global SU(2) algebra. 
Using the operators in Eq. (|117[) the Hamiltonian for a neutrino propagating through matter takes the form 



d 3 p- 



5m 



cos2(9J (p) + -sin26>(J+(p) + J_(p)) 



V2G F J d 3 pN e J (p). 



(118) 



In Eq. (|118l) , the hrst integral represents the neutrino mixing and the second integral represents the neutrino forward 
scattering off the background matter. Neutrino-neutrino interactions are described by the Hamiltonian 



y"d 3 pd 3 q(l-COstf pq ) J(p) • J(q) 



(119) 



where $ pq is the angle between neutrino momenta p and q and V is the normalization volume. Inclusion of antineu- 
trinos in Eqs. (|118j) and (|119[) introduces a second set of SU(2) algebras. For three flavors one needs two sets of SU(3) 
algebras, one for neutrinos and one for antineutrinos. Collective neutrino oscillations resulting from these equations 
exhibit a number of interesting symmetries [49 52] . 



5. CP-Violation in Neutrino Sector 



The neutrino mixing matrix is parameterized by three mixing angles and a CP-violating phase: 

/ 1 \ / C 13 S 13 e~ iS \ / C 12 Sia \ 
T 23 T 13 T 12 = (7a Sjg 1 -S 12 C 12 (120) 

V -S-23 C 23 J V S 13 e tS C13 / V 1 / 

where Cy = cos^-, Sij — sin^y, and S is the CP-violating phase. Only a non-zero value of #13 would also make the 
observation of the effects that depend on the CP-violating phase possible. Earlier hints for a non-zero value of #13 
from solar, atmospheric, and reactor data [53LI54I are further strengthened by the recent low-threshold analysis of the 
Sudbury Neutrino Observatory measurements [55| . Ongoing reactor experiments [5H4581 ] will provide a better insight 
into the value of this quantity. 

To explore the impact of neutrino propagation through matter on CP- violating effects we introduce the operators 

M 

* M = cos(9 23 1V - sin0 23 1 , T, 

* r = sin0 2 3*A.+ cos 23 * T , 
and write down the neutrino evolution equations as 

'!(•:)-=(»:) 

where 

(E x \ / Ve.fi \ 

H = T 13 T 12 E 2 T\ 2 T\ 3 + S$ 3 V Tli -C 23 S 23 V T ^ . (122) 
V E 3 J \ -C 23 S 23 V Tt , Cl 3 V T „ J 

In writing Eq. (I122p a term proportional to identity is dropped by adding a term to all the matter potentials so that 
= 0. Loop corrections in the Standard Model yield small, but non-zero values of V efl and V Tfl [60(. If we can 
neglect these terms it is straightforward to show that 

H(6) = SH(S = 0)S 1 

with 

S 

This factorization gives us interesting sum rules: Electron neutrino survival probability, P{v e — > v e ) is independent of 
the value of the CP-violating phase, 5; or equivalently the combination P(y^ — > u e ) + P(v T v e ) at a fixed energy is 
independent of the value of the CP-violating phase [6l| . It is possible to derive similar sum rules for other amplitudes 
(62j. These results hold even if the neutrino- neutrino interactions are included in the Hamiltonian [63j |. 
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